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Variational  Principles  for  Some^ 
Nonstandard  Elastic  Problems 

Variational  principles  are  derived  for  some  nonstandard  problems  involving  elastic 
bodies  in  smooth  contact.  For  these  problems,  the  portions  of  the  surfaces  where 
one  boundary  condition  holds  rather  than  another  must  be  determined  as  part  of  the 
solution  to  the  problem.  Cases  considered  include  a  body  containing  a  crack  or 
delamination,  indentation  bv  a  rigid  punch,  and  contact  with  an  elastic  foundation. 


Introduction 

The  principles  of  minimum  potential  energy  and  minimum 
complementary  energy  for  infinitesimal  elastic  deformations 
are  well  known  for  problems  in  which  the  tractions  or  the 
displacements  are  specified  at  each  point  of  the  surface  of  an 
elastic  body  (see  Sokolnikoff,  1956,  for  example).  The  prin¬ 
ciples  can  lead  to  bounds  on  quantities  of  physical  interest  and 
can  be  used  to  obtain  approximate  analytical  and  numerical 
solutions.  For  some  problems,  the  portions  of  the  surface 
where  one  boundary  condition  holds  rather  than  another  must 
be  determined  as  part  of  the  solution  to  the  problem.  For  ex¬ 
ample,  in  problems  involving  contact  between  elastic  bodies, 
the  shape  of  the  contact  area  can  vary  with  the  loading  (for 
references  see  Gladwell,  1980).  Uniqueness  of  solution  for 
typical  problems  of  this  type  was  considered  by  Shield  (1982) 
for  problems  involving  smooth  contact  between  surfaces  of 
elastic  bodies.  Flere  we  again  consider  elastic  problems  involv¬ 
ing  smooth  contact  and  develop  variational  principles  for 
some  typical  situations:  loading  of  a  body  containing  a  crack 
across  which  there  is  no  cohesion,  bodies  in  smooth  contact, 
indentation  by  a  smooth  rigid  punch,  and  contact  with  a 
smooth  elastic  foundation.  The  examples  can  be  combined  to 
treat  more  complex  problems,  such  as  the  indentation  of  an 
elastic  body  containing  a  crack. 

The  principles  rest  on  the  positive-definiteness  of  the  strain 
energy  and  they  show  that  the  potential  and  complementary 
energies  attain  absolute  minimum  values  only  when  the  trial 
functions  generate  the  strains  or  stresses  of  the  actual  solution. 
Weaker  stationary  principles  apply  without  the  assumption  of 
positive-definiteness  of  the  strain  energy.  It  is  assumed  that 
the  integrals  involved  are  convergent  if  they  are  improper  and 
this  requires  the  states  considered  to  have  finite  total  strain 
energy  or  total  complementary  energy.  The  elastic  material 
can  be  inhomogeneous.  At  an  interface  between  two  materials 
in  a  composite,  the  displacements  and  traction  are  assumed  to 
be  continuous  across  the  interface  but  delaminations  can  be 
included  if  they  are  modeled  as  cracks  across  which  there  is  no 
cohesion. 


Conmbuted  by  ihe  Applied  Mechanics  Division  for  publicaiion  in  ihc  Jour¬ 
nal  of  Applied  Mec  hanics. 

Discussion  on  this  paper  should  be  addressed  to  the  Editorial  Department, 
ASME.  United  Engineering  Center,  345  East  47th  Street,  New  York,  N  Y. 
10017,  and  will  be  accepted  until  two  months  after  final  publication  of  the  paper 
itself  in  the  Journal  of  Applied  Mechanics.  Manuscript  received  by  ASME 
Applied  Mechanics  Division,  January  2,  1987. 


Elastic  Body  With  a  Crack 

We  assume  that  the  strain-energy  density  M  of  the  body  is 
positive-definite  and  we  write 

2H'(e)=  c,jkle,jeu  ( c„k,  =  ckll/  =  c,lkl ) , 

where  etJ  are  the  infinitesimal  strains  referred  to  rectangular 
Cartesian  axes  x,.  The  stresses  ttJ  are  related  to  the  strains 
through 

A  =  C ijkl^kl ’  Gy  ~  Cy^./Cy, 

where  C,,kl  have  the  same  symmetries  as  c,/A/.  The  strain- 
energy  is  also  a  positive-definite  function  W(  of  the  stresses, 

2  Wc  (t)  =  Cljklt,jtkl. 

In  the  unstressed  reference  state,  the  body  occupies  a  region 
V  with  surface  S  and  we  suppose  that  it  contains  a  crack  across 
which  the  material  has  no  cohesion.  The  crack  is  defined  by  a 
surface  C  in  K,  and  we  use  n  to  denote  the  unit  normal  to  one 
side  of  the  surface  C.  We  use  ±  signs  to  indicate  values  of 
quantities  on  the  two  sides  of  C,  with  the  +  sign  referring  to 
the  side  of  C  with  exterior  normal  n.  We  shall  also  use  square 
brackets  to  indicate  the  difference  in  the  values  of  a  quantity 
across  C,  so  that  for  the  displacement  field  u,  for  example, 

[u]  =  u  *  -  u  . 

Under  loading  of  the  body,  we  assume  that  at  points  of  the 
crack  surface  C  either  (/)  the  crack  opens  with  no  traction 
transmitted  across  C  or  (it)  the  two  sides  remain  in  smooth 
contact.  Then  at  points  of  C  we  have 

(u  ‘  -  u  )*n  =  [u]*n<  0,  •  (1) 

and  we  require  on  C 

either  (i)  T *  =  —  T  =0  when  [u]-n  <  0  T 

.'A  (2) 

or  (//)  T*  =  -  Tf  =  -pn  when  fu]  =  0,  J 

b 

where  T  is  the  surface  traction  and  p(x)  is  the  (nonnegative) 
pressure  transmitted  across  C.  We  set  p-  0  at  points  of  C 
where  (/)  holds. 

The  loading  of  the  body  is  assumed  to  be  caused  by  a  body 
force  F.  prescribed  surface  tractions  T''  on  a  portion  Sr  of  S 
and  prescribed  displacements  u<;  on  the  remainder  S(  of  S. 
The  displacement  field  u  is  unique  except  possibly  for  a  rigid 
displacement,  depending  on  the  conditions  on  S,  (Shield, 
1982). 

We  define  the  potential  energy  P  as  the  functional 
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P\U'\  = 


\  W(e’)d V-  \  T<;*u'c/.S-  i  F-uVF 

J  t  Js,  Jl 


(3) 


for  fields  u '  such  that 

u '  =  uG  on  SL  ,  [u ' ] *n  <  0  on  C. 

We  set 


u'  =  u  +  Au. 


(4) 


Then  Au  is  zero  on  SL  and  we  have 

H'(e')=  W{e)  +  t„Ae„  +  W'(Ae).  (5) 

The  stresses  itJ  satisfy  equilibrium  with  body  force  F  and  with 
the  divergence  theorem  and  boundary  conditions  we  get 


P|u'|-P(w|  = 


W'(A e)dV+  \  T*  [Au]rfS, 
i  .'c 


using  the  divergence  theorem  and  T  •  [u]  =  0  on  the  crack 
(Clapeyron’s  theorem).  We  then  have 

Pfw'l  >  P[u\  >  -Q[f  |,  (9) 

and  the  principles  provide  upper  and  lower  bounds  for  the 
potential  energy. 

Smooth  Contact  Between  Elastic  Bodies 

For  simplicity  we  consider  contact  between  two  elastic 
bodies  occurring  over  parts  of  the  surfaces  of  the  bodies  which 
are  nearly  plane,  but  the  approach  is  easily  generalized.  In  the 
reference  configuration  the  bodies  touch  at  the  origin  0  of 
coordinates  and  the  plane  x}  =  0  is  tangent  to  both  bodies  at 
0.  The  bodies  occupy  regions  V,  and  Vz  with  the  a,  axis 
pointing  into  K:.  Under  loading,  contact  may  occur  over  sur¬ 
faces  C |  and  C,  of  the  bodies,  defined  as  the  nearly  plane  sur¬ 
faces 


where  contributions  from  both  sides  of  C  have  been  included 
in  the  integral  over  C.  For  the  actual  solution  T  •  [u]  is  zero  on 
C  in  view  of  equations  (2)  and  the  integral  over  C  has  the  value 

^  T*[u']tfS=  -  ^  pn*(u'|rfS. 

r;om  equation  (4)  this  is  seen  to  be  greater  than  or  equal  to 
zero  and  with  IV  positive-definite,  we  have 

P(u'  |  >  P\u\ 

with  equality  if  and  only  if  u  and  u  '  have  the  same  strains  and 
n  •  [u‘]  is  zero  where  p  is  nonzero. 

Thus  we  have:  For  displacements  which  satisfy  the  displace¬ 
ment  boundary  conditions  and  have  no  interpenetration  of 
material  across  the  crack ,  the  potential  energy  P  is  least  for  the 
displacements  of  the  actual  solution. 

The  complementary  energy  Q  is  defined  as 

Q\f t  =  Wc(t')dV-  T'  •ucdS  (6) 

for  stresses  in  equilibrium  with  F  and  such  that 

=  Tf  on  Sr,  T'  *  =  -T' '  =  - p'n  on  C, 
where p'  >  0.  We  set 

so  that  the  stresses  A r,y  satisfy  equilibrium  with  no  body  force 
and  have  zero  traction  on  ST.  We  have 

Wr  ( t ' )  =  Wc  ( t)  +  A  t,;e„  +  fT(Ar),  (7) 

and  with  the  divergence  theorem  we  obtain 

Q\t'  I  ~QU\  =  j(  Wc  (A l)dV+  Ar*(ulrf5. 

As  before  T  •  [u ]  is  zero  on  C  and  the  integral  over  C  becomes 

j  T'»[u]cfS=-  J  p'n*[uldS  >  0. 

It  follows  that 

Q\f  I  2  Q\t\ 

with  equality  if  and  only  if  t,'  =  ttJ  in  V,  and  we  have:  For 
stress  fields  with  the  given  surface  tractions  and  in  equilibrium 
with  the  given  body  force  and  which  transmit  at  most  pressure 
across  the  crack,  the  complementary  energy  Q  is  least  for  the 
stresses  of  the  actual  solution. 

For  the  actual  solution  we  have 

P\u\  i-Q!/)  =2  f  WdV-  \  T'udS-  {  F*uc/F  =  0,  (8) 


C,:  .v3  =/(.*,,  .v2).  C::  .v,  =g(.v,,  .v;), 

where  a  , ,  x2  lie  in  a  region  C  of  the  .v,  -  as  plane  enclosing  the 
origin  and 

/(at,,  a2)  <  g(A-,,  a3). 

The  contact  is  smooth  and  we  denote  the  pressure  between 
the  two  bodies  by  the  nonnegative  function  p(x, ,  ,v;),  defined 
over  C.  Then 


p\u;2 -u;1  -  (f-g )  i ds, 

and  this  is  nonnegative  because  uj  satisfies  (II).  Thus  we  can 
show  that:  For  displacements  which  satisfy  the  displacement 
boundary  conditions  and  which  satisfy  (II)  on  C,  the  poten¬ 
tial  energy  is  least  for  the  displacements  of  the  actual  solution. 
The  complementary  energy  Q  is  defined  by 

=  (  W,.(r)dV-\  T'.u«rfS-j  p'  (f-g)dS  (13) 

•I  *  J  S(j  Jc 

for  stresses  r ^  in  equilibrium  with  the  given  body  force  and 
given  tractions  and  which  involve  at  most  a  pressure  p'  be¬ 
tween  the  bodies  across  C.  The  difference  between  Q\t' )  and 
01  f !  is  transformed  as  in  the  previous  section  and  we  are  led 
to  consider  the  sign  of 


T\  =  -  T\  =p(*|,  .v:)  on  C,  (10) 

where  the  superscripts  indicate  values  for  the  two  bodies.  The 
bodies  do  not  penetrate  each/so  that  the  displacement  compo-  .  ,, 
nent  m,  satisfies  ^  ^ 

u]-u\  >  f-g  on  C.  (11) 

At  each  point  of  C  we  require 

p  =  0  when  -  u\  >  f-g ,  u]-u\=f—g  when  p  >  0.  (12) 

On  the  remaining  portions  5,  and  S2  of  the  surfaces  of  the  two 
bodies,  we  suppose  that  tractions  TG  are  specified  on  parts 
SlT,  S2T  and  displacements  uc  are  specified  on  parts  S,t ,  S:L, 
of  S, ,  S;,  respectively.  We  use  Sr  and  SL  to  denote  S17  +  S;r 
and  S,(;  +  Szu.  The  body  force  F  is  assumed  known  in 
V  =  V,  +  Vz.  The  solution  is  unique  except  possibly  for  a  rigid 
body  displacement  (Shield,  1982). 

The  potential  energy  P  is  again  defined  by  equation  (3)  for 
displacements  which  satisfy  the  boundary  conditions  on  Sc 
and  which  satisfy  (II)  on  C.  As  in  the  previous  section,  we 
look  at  the  difference  between  P\u'  \  and  P\u\,  use  equation 
(5)  and  the  divergence  theorem,  and  find  that  the  difference 
involves  the  integral 

\  T*(u '  -  u)dS. 

JCftCy 

In  view  of  equations  ( 10)  and  (12),  the  integral  becomes 
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(18) 


(T'-T).udS-  (/?'  -p){f-g)dS. 

,  *  c\  Jr 


P\u'  |  =  ^  W(e')dV-  T(- 


•  u '  dS 


This  can  be  written  as 

f  (P'  ~P)  (tn -«S -W~g) \dS. 

J  c 

and  with  equations  (11)  and  (12)  we  see  that  the  integral  is 
nonnegative.  Thus  we  can  show  that:  For  stress  fields  in 
equilibrium  with  the  given  body  force  and  given  tractions  and 
wh’ch  involve  at  most  a  pressure  between  the  bodies  across  C, 
the  complementary  energy  is  least  for  the  stresses  of  the  actual 
solution. 

For  the  actual  solution  we  again  have  P\u\  =  -  Q\  t )  and 
the  principles  lead  to  upper  and  lower  bounds  on  P|tr|. 


-  F*u  dV-Ld'  -A/,0' 

for  fields  u'  which  satisfy  the  displacement  conditions  on  S0 
and  which  are  such  that 

u'3  <  g(.V|,  x2)-d ' +o '.v:  -  b '.v,  on  C,  (19) 

where  d' ,  o' ,  b'  are  constants.  After  transformation,  we  find 
that 

P|u' ) -P\u\  =  l  e)dV 


Indentation  by  a  Smooth  Rigid  Punch 

In  order  to  illustrate  problems  in  which  an  elastic  body  can 
come  into  contact  with  a  rigid  body  of  known  shape,  we  con¬ 
sider  indentation  by  a  smooth  rigid  punch  when  the  possible 
area  of  contact  is  a  region  C  of  the  x,-*.  plane  enclosing  the 
origin  0.  The  exterior  normal  to  the  body  at  0  is  along  the  x, 
axis,  and  the  remainder  of  the  surface  of  the  body  is  denoted 
by  S. 

When  the  movement  of  the  punch  is  known,  we  will  have 

n,  <  g(x,,  x2)  on  C,  (14) 

where  g  is  a  known  function,  and  we  require  at  each  point  of  C 
p  =  0  when  n,  <  g,  u}=  g  when  p  >  0,  (15) 

where  p  is  the  contact  pressure.  Tractions  are  prescribed  on  a 
part  ST  of  S  and  displacements  on  the  remainder  SL  of  S,  with 
a  known  body  force  in  the  region  V  occupied  by  the  body. 

The  problem  can  be  considered  as  a  limiting  case  of  contact 
between  two  elastic  bodies.  The  potential  energy  is  defined  to 
be  the  functional  (3)  for  displacements  which  satisfy  the  condi¬ 
tions  on  St.  and  (14)  on  C.  We  can  then  proceed  as  in  the 
previous  section  to  show  that:  The  potential  energy  is  least  for 
the  displacements  of  the  actual  solution. 

The  complementary  energy  Q  is  defined  to  be 

(?(?')  =  Wc(t')dV-  T'.uc;c(S-f  p'gdS 

for  stresses  in  equilibrium  with  the  given  body  force  and  sur¬ 
face  tractions  and  which  involve  a  pressure  p'  in  the  contact 
area  C.  Then:  Q  will  be  an  absolute  minimum  for  the  stresses 
of  the  actual  solution. 

Instead  of  prescribing  the  movement  of  the  punch,  we  may 
prescribe  the  downward  force  L  on  the  punch  and  the 
moments  A/,,  M2  of  the  force  about  the  x,,  x2  axes,  with 
prescribed  loading  on  S  as  before.  The  contact  pressure  p  must 
then  satisfy 

L  =  ^  pdS,  A/,  =  -  j  px2dS ,  A/2  =  j  Px\ dS.  (16) 

For  a  known  punch  shape  g(x, ,  x2)  the  solution  will  satisfy 

Uj  <  g(.V|,  x2)-d  +  ax2 -bx,  on  C,  (17) 

with  equality  where  p  is  nonzero.  The  constants  d,  a ,  b  are 
determined  as  part  of  the  solution,  which  will  be  unique  except 
possibly  for  a  rigid  body  displacement  (Shield,  1982).  (Other 
problems  may  be  treated;  for  example,  we  may  require  the 
punch  to  indent  without  tilting  and  then  a,  b  are  zero  and  A/, , 
M2  are  not  prescribed.) 

The  potential  energy  is  defined  as 


p(u(  - u2)dS—  L{d'  -d)  -Mfa'  -a)  -M2{b‘  -b). 

If  we  set 

v  =  u j  +d  —  ax2  +  bxt,  v'  =  u{  +  d’  —a'x2  -t -b'xx 

and  use  equations  (16),  then  apart  from  the  strain-energy  term 
the  right-hand  side  becomes 

p(v’  —  v)dS. 

Now  where  p  >  0,  equality  holds  in  (17)  and  v  =  g.  Because 
v'  <  g  from  (19),  we  then  see  that  the  integrand  is  non¬ 
positive.  Thus:  The  potential  energy  is  least  for  the 
displacements  of  the  actual  solution. 

The  complementary  energy  is  defined  to  be 

Q[f\=\  Wc(t')dV-[  T'»ucdS+  (  p’gdS 

d  V  J  S(j  J  c 

for  stresses  in  equilibrium  with  the  given  body  force  and  sur¬ 
face  tractions  and  which  involve  a  pressure  p‘  in  the  contact 
area  C  satisfying  the  loading  conditions  (16).  We  can  then 
show  that 

Q\T  1  ~Q\t\  =  Wc(At)dV-  j  (p’  - p){u ,  -g)dS. 

Because  p'  and  p  apply  the  same  resultant  force  and  moment 
over  C,  the  integral  over  C  can  be  written  as 

j  {p‘  —p)  (u3  -  g  +  d—ax2  +  bxfdS. 

Equality  holds  in  (17)  wherep  >  Oso  that  there  is  no  contribu¬ 
tion  to  the  integral  from  p,  and  the  contribution  from  p'  is 
nonpositive  in  view  of  (17).  It  follows  that:  The  complemen¬ 
tary  energy  is  least  for  the  stresses  of  the  actual  solution. 

For  the  actual  solution  we  have  P|w)  =  —  Q\t\  for  both 
punch  problems  of  this  section. 


Smooth  Contact  With  an  Elastic  Foundation 

Variational  principles  also  hold  for  elastic  bodies  which  can 
receive  support  from  a  foundation  of  the  Winkler  type.  For 
simplicity  we  assume  here  that  the  foundation  has  a  plane  sur¬ 
face  and  lies  below  the  plane  x,  =  0.  The  petition  of  the  sur¬ 
face  of  an  elastic  body  that  can  come  into  contact  with  the 
foundation  is  the  nearly  plane  surface 

x\  =f(x,.x2) 

touching  the  plane  x2  =  0  but  lying  entirely  above  it,  where 
x,,jc2  lie  in  a  region  C  of  the  .x^.v.  plane.  The  reactive 
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pressure  p  of  the  foundation  is  proportional  to  the  downward 
displacement  of  the  surface,  so  that  on  C  we  require 

p  =  0whenu!  >  —f,p- -K(u^  +f)  when  u,  <  -/,  (20) 

where  A.'  >  0  is  the  stiffness  of  the  foundation  (A  may  vary 
with  x, ,  x:).  if  we  define  q(u,)  by 

q(ui)  =  0  when  u}  s  -/,  q(u})  =  1  when  u,  <  — /, 
then  we  can  write 

7",  =  /?  =  -K(u}  +f)  q(U))  on  C,  (21) 

while  the  tangential  tractions  are  zero  on  C.  Boundary  condi¬ 
tions  on  the  remaining  surface  S  of  the  body  and  a  body  force 
field  are  prescribed  as  before. 

The  potential  energy  is  defined  to  be 

P\u  )  =  ^  W(e')dV-  ^  T°*u  dS 

-f  F.UW+—  [  K(u{+f)2q'dS, 

J  t  2  J  c 

where  q  '  =  q(uj),  for  displacement  fields  //,'  which  satisfy  the 
displacement  boundary  conditions  on  St,.  Using  equation  (S) 
and  the  divergence  theorem,  we  find  that 

Pin'  |  -P\u\  =  1  W(Ae)dV 

J  V 

+  -y  K(v'2q'  +  v2q-2vqv' )dS  (22) 

in  which 

V  =  U}+f,  v'=U)+f. 

By  considering  the  various  possibilities  for  the  signs  of  v  and 
v‘  together  with  the  corresponding  values  for  the  step  func¬ 
tions  q  and  q’ ,  we  find  that  the  integrand  of  the  integral  over 
C  in  equation  (22)  is  nonnegative.  Thus  we  again  have:  The 
potential  energy  is  an  absolute  minimum  for  the  displacements 
of  the  actual  solution. 


For  the  complementary  energy  we  take 
Q\C  I  =  Wc(f)dV-  ^  T'-u<;rfS 


+  J  [p'2/2K+p'f)dS 

for  stresses  in  equilibrium  with  the  given  body  force  and  sur¬ 
face  tractions  and  which  involve  a  pressure  p'  on  the  founda¬ 
tion  inte-face  C.  Using  equation  (7)  and  the  divergence 
theorem,  we  can  show  that 

■*.i't-QWI  =  Wc{\t)d\/ 


+  jt.  I  (p'2 -pz)/2K  +  [p'  - p){u ,  +f) \dS. 

From  equations  (20),  the  integrand  of  the  integral  over  C  is 
found  to  be  nonnegative  and  we  see  that:  The  complementary 
energy  is  an  absolute  minimum  for  the  stresses  of  the  actual 
solution. 

We  can  also  show  that  P|w|  =  -  Q\t\,  so  that  the  prin¬ 
ciples  can  be  used  to  bound  P(  u  | . 
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